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A SIMPLIFIED METHOD FOR PREDICTING THE STABILITY OF 
AERODYNAMICALLY EXCITED TURBOMACHINERY 

Albert F. Storace 
General Electric Aircraft Engines 
Cincinnati, Ohio, U.S.A. 


A method is presented for the quick and accurate prediction of the stability of 
aerodynamically excited turbomachinery using real eigenvalue/eigenvector data 
obtained from a rotordynamics model. An expression is presented which uses the 
modal data and the transmitted torque to provide a numerical value of the 
relative stability of the system. This approach provides a powerful design tool 
to quickly ascertain the effects of squeeze-film damper bearings, bearing 
location, and support changes on system stability. 

INTRODUCTION 

The purpose of this paper is to present a method that is easily and economically 
applied to turbomachines to predict the effects of shaft flexibility, 
squeeze-film bearing supports, and static structure configuration on the 
rotor-bearing/static structure system stability relating to rotor aerodynamic 
cross coupled stiffness (Alford) forces. The method is general in that systems 
with general rotor support arrangements and multiple spools can be handled. A 
major advantage of the method is that it allows machine designers to quickly 
determine the effects on stability of bearing changes, shaft modifications, and 
bearing support designs to determine appropriate system designs. This paper 
presents an expression and analysis methodology for predicting system stability 
that includes the effects of destabilizing forces, rotor/stator dynamic 
displacements, internal and external damping, and gyroscopic moments. The 
expression and methodology presented provides an analysis approach that is 
simplified but at the same time includes all of the modeling detail needed to 
perform a valid assessment of system stability. 

The method provides a timely and cost effective means to initially screen designs 
without incurring the high computer costs and large amounts of data reduction 
time required using complex eigenvalue rotor dynamic analysis programs in a 
repetitive mode. Specifically, the method uses the results of an undamped 
lateral critical speed analysis and generalized forces derived from the physical 
destabilizing forces to develop modal equations of motion for a self-excited 
system. The solution of these equations is then used to develop a modal 
dimensionless stability criterion. This criterion requires that the energy 
absorbed by the system exceed the energy imparted to the rotor by unbalanced 
torque forces if instability is to be prevented. 

The method is intended to supplement the more general analysis techniques which 
are used for final confirmation of the stability predictions. 
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SYMBOLS 


F X’ F XR’ F Y’ f yr 


F XS’ F YS 


X, Y, X R , Y r 

X S’ Y s 
K XY’ k yx 
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P X’ P Y 


^XR’ ^YR’ ^XS’ ^YS 


G X’ g y 
m 

C 


il) 

[K] 

CO 

Q 

N 

RPM 

nb 


Aerodynamic forces acting on the rotor N 
(LBf) 

Aerodynamic forces acting on the stator N 
(LBf) 

Rotor displacements cm (IN.) 

Stator displacements cm (IN.) 

Cross-coupling spring rates N/cm (LBf/IN.) 

Compressor or turbine stage torque N-cm 
(LBf-LB) 

"Alford" coefficient (dim.) 

Stage pitch diameter cm (IN.) 

Blade height cm (IN.) 

Generalized coordinate pair cm (IN.) 

Modal displacements (dim.) 

Generalized forces N (LBf) 

Generalized mass Kg (LBf-SEC 2 /IN. ) 

Generalized damping coefficient N-SEC/cm 
(LBf-SEC/IN. ) 

Generalized stiffness coefficient N/cm 
(LBf/IN.) 

Mode shape vector (dim.) 

Physical stiffness matrix N/cm (LBf/IN.) 

Undamped natural frequency (RAD/SEC) 

Modal Q-factor (dim.) 

Number of rotor stages (dim.) 

Design point speed (REV/MIN) 

Number of field damping components 
(casing, frames, rotors) 
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SYMBOLS (Cont.) 


na Number of lumped damping components 

(mounts, bearings, dampers) 

Qa Component Q-factor for lumped damping 

components (dim.) 


Qb 




N 


1 


N 


2 

MSC 


Component Q-factor for field damping 
components (dim.) 

Spin speeds for field damping components 
(#. is zero for static components) 
(RAO/SEC) 

Low pressure rotor speed (REV/MIN) 

High pressure rotor speed (REV/MIN) 

Modal Stability Criterion (dim.) 


BACKGROUND 


A major destabilizing mechanism acting on turbomachinery stages is the Alford 
aerodynamic cross-coupling stiffness force (ref. 1). In a fixed frame global 
coordinate system, this force can be modeled by the following equation. 
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Figure 1, the 


( 1 ) 


(T/3/DpH) N/cm, T is the stage torque, D p is the 
of the '^tage','' H is the Dlade height, and /3 is the change in 
efficiency per unit change in blade tip clearance, expressed as 
blade height. The physical rationale for these forces is based on 
blade efficiency with decreasing tip clearance. Referring to 
displacement of the disc centerline resulting from whirl decreases the blade tip 
clearance in the direction of the displacement. The efficiency of those blades 
with reduced clearance is improved resulting in a greater than average torque 
delivered by those blades with reduced clearances. Conversely, on the side of 
the disc with increased clearances, a less than average torque is imparted to the 
rotor by those blades. The integrated effect of the circumferential variation of 
blade torque results in a net torque in the direction of whirl associated with 
the vector force described in equation 1. As discussed in references 1 and 2, 
the cross coupled stiffness can be responsible for self-excited rotor instability 
at high power levels that is characterized by subsynchronous rotor whirl in the 
direction of rotation. This subsynchronous whirl is generally associated with the 
first rotor dominated mode and the stability analysis method to be presented in 
this paper addresses a mode by mode evaluation approach to ascertain the 
potential for system instability. 
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STABILITY CRITERION 


Equation 1 expresses the physical forces acting on the rotor. It can be extended 
to include the forces acting on the case (stator) as follows. 
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These forces will be used in conjunction with the gyroscopical ly stiffened modes 
obtained from a real eigenvalue/eigenvector analysis to develop a stability 
criterion. For a given mode obtained from the real mode set, define Pw and Pw as 
the modal coordinate pair describing the generalized response in the vertical 'and 
horizontal planes of the system. More will be said later in this paper about the 
incorporation of gyroscopic effects in the modal data. 


The rotor and stator modal displacements at a given stage can then be used to 
define the physical destabilizing forces acting on the rotor and stator as 
follows: 
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(3) 


The total generalized forces acting in the X and Y directions can be written 


G x " " k xy p y ^yr'M ^xr + k xy p y ^yr " M hs _ 

= -K XY P y (0 yR - 0 ys ) (0 XR - 0 X$ ) 

G Y = K YX P X ( ?XR ' ^YR ' K YX P X ( ^XR " W ^YS = 
= K yx P x (0 xr - 0 XS ) (0 yR - p ys ) 

If circular whirl is assumed, 



(5) 


Then the generalized forces are 



( 6 ) 


The equations of motion for the modal coordinate pair are 
mPy + CPy + k Py = Gy 

m*P Y + CP^ + k P(J = G} (7) 
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Combining these equations and eliminating Py, 
d^P x 2C_ d 3 P x /2k C 2 ^d Z P x 2Ck dPj 

3 +( ' 


dt 


m dt 
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.st 


For a solution, assume P x = p e"”' and upon substitution, the following 
characteristic equation is obtained. 


S 4 + A 3 S 3 + A 2 S 2 + AjS + A q 


( 9 ) 


The complex eigenvalues of this equation define the stability boundaries of the 
system. Specifically, a root with a positive real part indicates an unstable 
system. Per the Routh stability analysis, an inspection of the coefficients of 
the characteristic equation determines whether the motion is stable or unstable 


For a stable system, 


A i A 2 A 3 > A | + Ag Aq 


(10) 


Evaluating the terms, 
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Equation 10 can be written as 
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where CO =^|k/ m = undamped natural frequency, 
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and Q = k/COO = modal Q-factor. 

Expressing the relative modal displacement on a per stage basis 


and the cross-coupling stiffness in terms of the rotor speed and 
the HP per stage at the design point leads to the following equation. 


For stability. 


2PE 


713,361 N A. HP, 

q — — — (* R1 

RPM 1-1 D pi H i K1 



> 1.0 ( 12 ) 


T 

where PE = 1/2 {0} [K](0> = 1/2 k = modal potential energy and N = the number of 
rotor stages. Equation 12 represents an energy balance expressed in terms of the 
pertinent modal parameters for a system mode of vibration and the physical 
destabilizing forces. If English units are used, then the constant 713,361 
RPM-cm-N/HP becomes 63,025 RPM-IN-LB/HP. The accurate calculation of the modal 
or generalized Q-factor is key to the use of equation 12 for evaluating system 
stability. It must reflect the effects of both external and internal damping and 
the modal participation of the various engine components. 


The modal Q-factor is given by equation A10 in the Appendix. Substituting this 
equation into equation 12 yields the following modal stability criterion (MSC). 
For stability, 
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> 1.0 (13) 
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Note that squeeze-film damper elements contribute lumped damping and are included 

in the summation a-1, 2, 3, na. As an approximation, they are modeled as 

soft springs in the system vibration analysis and a conservative component 
Q-factor of 3.0 can be used, although a more exact value can be calculated. 


GYROSCOPIC EFFECTS AND ROTOR INTERNAL DAMPING 


Figure 2 shows an example of a typical engine system vibration model used to 
generate modal data for the MSC. This model represents a single plane of a 
demonstrator engine and is an assemblage of substructure (span) and spring-type 
elements. It can be built very rapidly and is easily altered and interactively 
run to generate modal data for MSC evaluation of a wide range of alternative 
system designs. The span element type, represented by solid lines in Figure 2, 
includes both flexibility and mass properties and models casings, rotors, and 
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frames. The spring-type elements model bearings, mounts, and dampers. 
Gyroscopic moments are incorporated through spin and whirl frequency dependent 
terms in the mass matrices of the substructures. Note that while the modal data 
is obtained from an analysis which models a single plane of the engine, it does 
reflect whirling motion of the rotors. Consider the consequence of cross-axis 
gyroscopic coupling: (1) lateral motions of the rotor are not planar-the rotor 
center motion describes a circular orbit, if rotational symmetry prevails; (2) 
each free vibration mode of the equivalent non-rotating shaft of the planar model 
is split into two modes which are distinguishable by the sense of whirl motion 
(relative to the shaft spin). These are forward and backward whirling modes, and 
since the Alford instability mechanism drives forward whirling modes, the planar 
model is constrained to provide forward whirling modes for the reference rotor. 
Figure 3 shows the Campbell diagram (map of natural frequencies vs spin speeds), 
for the model of Figure 2, referenced to high pressure rotor spin speeds. The 
frequency lines represent system modes involving forward whirl of the high 
pressure rotor and backward whirl of the low pressure rotor. The latter are a 
consequence of counterrotating rotors, and decreasing natural frequency (due to 
gyro softening) with increasing high pressure rotor speed reflects dominant low 
pressure rotor participation. Figure 4 shows examples of mode shapes for two 
high pressure rotor subsynchronous modes at the 13,226 RPM N«/ - 1 1 ,340 RPM N, 
design point for the Campbell diagram of Figure 3. The SlSl^cycle/min (CPM) 
mode is a fan shaft bending mode and the 8303 CPM mode is a core rotor bending 
mode, with the core rotor out-of-phase with. the core case. For these two modes, 
the spin speed to whirl frequency ratios jfl./co f° r the high and low pressure 
rotors are equal to (4.22, -3.62) and (1.59, -1.36), respectively. The MSC 
evaluations are based on a mode by mode evaluation at the design point (design 
speeds and reference rotor torque). The MSC values are calculated for each 
forward whirl high pressure rotor mode and the minimum value is the basis for the 
rotor system stability prediction. Since the destabilizing effects of the 
aerodynamic cross-coupling forces are generally much more significant for the 
high pressure rotor than for the low pressure rotor, the former is considered as 
the reference rotor. Hence, the index i in equation 13 ranges over the stages 
for the high pressure rotor; both power absorption (compressor stages), and power 
generation (turbine stages) are included in the summation. It will be noted that 
single mode evaluation is acceptable because the cross-axis stiffness associated 
with the Alford forces is relatively small and, therefore, little loss in 
accuracy results from the use of the original mode shapes in the stability 
calculations. 

The incorporation of the modal Q-factor in the MSC results in the implicit 
inclusion of rotor hysteresis or rotary damping which can be destabilizing if the 
rotor is undergoing subsynchronous vibration. The spin speed to whirl frequency 
ratios obtained from the system vibration analysis provide the data needed to 
correctly incorporate the effects of rotor internal damping in the generalized or 
modal damping for the mode. Hence, the modal Q-factor provides the effective 
system damping needed to correctly define the energy absorbed by the turbomachine 
at resonance. 


EXAMPLE STABILITY CALCULATION AND ROTOR WHIRL EXPERIENCE 


To demonstrate the calculation method, the MSC values are calculated for each 
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natural frequency involving high pressure rotor forward whirl at the design point 

for the engine model of Figure 2. These natural frequencies are calculated up to 
the HP rotor synchronous frequency and correspond to the intersection of the 
Campbell diagram frequency lines and a vertical line passing through the 
reference rotor design speed (Figure 3). They occur at 1118, 1554, 1665, 3131, 
4034, 5966, 6654, 8303, 10632, and 12580 CPM. 


For example, for the 8303 CPM high pressure rotor bending mode, the numerator of 
equation 13 is equal to 122,708. In general, 1<<3<2. Setting (-3 equal to 2.0 
for each stage, the term 


1L @ i HP i 


i=l DP i H i 


^Ri ' /*Si) 


9 

of equation 13 is equal to 1569 HP/cm . 


As previously mentioned, the summation i = 1 , 2, N encompasses both the 

compressor and turbine stages of the high pressure rotor. Then at the = 
13,226 RPM design speed of the high pressure rotor: L 


The MSC is equal to 123,067 N/cm 

, = i .45 

713,361 

(1569) N/cm 

13226 


Table 1 provides the MSC and modal Q-factor for the design point modes. This 
table shows that the minimum value for the MSC occurs for the 8303 CPM mode and 
that the system is predicted to be stable. 

The modal stability criteria has shown good correlation with experience for 
various General Electric Aircraft engines. 


CONCLUSIONS 


The MSC provides a convenient and quick means to perform a rotor stability 
analysis using modal data readily available from planar system vibration models. 
It includes all of the significant parameters (gyroscopic moments, damping, 
rotor/stator relative displacements), and modeling detail needed to perform a 
valid assessment of rotor stability related to Alford forces. The assumption of 
circular whirl means that the stabilizing effects of non-axi symmetric rotor 
and/or engine support stiffness are not included, and this may result in built-in 
conservatism, depending on the characteristics of the engine system modes. 
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TABLE I - MODAL STABILITY AND Q-FACTOR VALUES CALCULATED AT THE 
DESIGN POINT* FOR THE ENGINE MODEL OF FIGURE 2 


NATURAL FREQUENCY 
(CPM) 

MODAL-Q 

MSC 

1118 

12.5 

1657 

1544 

w 

to 

2643 

1665 

15.1 

1112 

3131 

20.7 

12 

4034 

16.1 

58.9 

5966 

21.4 

246 

6654 

15.3 

7.5 

8303 

19.3 

1.45 

10632 

16.5 

53.7 

12580 

18.5 

6.0 


* N 1 = 11340 RPM/N 2 = 13226 RPM 
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CLEARANCE 



ROTOR DISK ECCENTRIC TO STATOR IN X DIRECTION 
PRODUCES NET FORCE IN Y DIRECTION 


FIGURE 1 - ALFORD INSTABILITY MECHANISM 



FIGURE 2 - SYSTEM VIBRATION MODEL OF A DEMONSTRATOR ENGINE 
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FIGURE 3 - CAMPBELL DIAGRAM 

FORWARD WHIRL MODES FOR HP ROTOR 
BACKWARD WHIRL MODES FOR LP ROTOR 
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HP ROTOR FORWARD WHIRLING MODE AT 
THE DESIGN POINT 


Nl=-11340, N2=13226 RPM=3131.3 

TOTAL ENERGY 128744 



R =13226=+4.22 
3131 


HP ROTOR FORWARD WHIRLING MODE 
AT THE DESIGN POINT 


Nl=-11340, N2=13226 RPM=8302.6 

TOTAL ENERGY 679347 



R =-11340= -1.36 
XP 8303 

r hp = 13226=4-1.59 
8303 

FIGURE 4 - EXAMPLES OF HIGH PRESSURE ROTOR SUB SYNCHRONOUS MODES 
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APPENDIX 


DEVELOPMENT OF MODAL Q-FACTOR EQUATION 


The equation for the modal Q-factor is developed by summing the damping 
contributions of the engine components (field damping for casings, frames, and 
rotor component structures, and lumped damping for mounts, bearings, and 
squeeze-film dampers). 

The Q-factor for the i-th component is defined as 

27TPE< 

Q h - (Al) 


where PE| is the physical strain energy in the component and Ep. is the energy 
dissipated. 


Consider the work done on the system at resonance by an excitation force P at 
point j expressed in terms of modal data for the system. 

"in ■ TP /j (SF) (fl2) 

where #. is the modal displacement at point j and SF is a scale factor relating 
the modal displacement to the physical displacement. 

From equation Al, the energy dissipated by the engine components 


is 

9 n PE. 

W n = 2 7T(SFr 2Z — ~ (A3) 

U i=l Qi 

where n is the number of components and PE. is the modal strain energy in the 
i-th component. 1 


Equating the work done on the system to the energy dissipated leads to a solution 
for the modal scale factor SF. 
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(A4) 
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where f- is the modal displacement at point j. 

J 

In reality, the modal scale factor is the modal participation at resonance, or 


SF = 


= 3 l 


(A5) 
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e*^ , 

where j: = P.y. = the generalized force and k is the generalized stiffness. 

J J 


Hence, k 

Q = (SF)^r 
cf 


P.0. 

J P J 


PE, 
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221 
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(A6) 


Equation A6 provides the modal Q-factor for the system reflecting the damping 
contributions of the static structures and the rotor component structures. 
However, in the latter case, the rotors are treated as stationary component 
structures since the effects of spin and whirl have not been included. Hence, 
the internal or hysteretic damping associated with a spinning and whirling rotor 
has been neglected. This damping mechanism is characterized by an internal 
friction force caused by the rate of change of strain within the rotor. This 
internal friction force can be represented with the viscous damping model derived 
in references 3 and 4 as 


F t = -C t cf(a>-J h (A7) 

where C. is the viscous damping coefficient for a stationary rotor, jS is the spin 
speed, ^>is the whirl frequency, and is the whirl displacement. 


Equation A7 represents the follower force (tangential force which leads the whirl 
displacement cT by 90°) F. in a rotating coordinate system fixed to the rotor. An 
equivalent damping coefficient is derived as follows. 

F t = -C eq £JcT = -C = -C^CQqT (l-jt/co) (A8) 

Hence, C £ q = C t (1 -$/a>) 

If Q r is the component Q-factor for a stationary rotor corresponding to C., then 
an equivalent field Q-factor Q™ for a whirling and spinning rotor of stiffness 
k^ can be defined as follows. 4 
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EQ = 


or 


1 _L 

q eq co q r co 

Q .-j£_ 

EQ (1-0/W) 


(1 -jl/a» 


(A9) 


Notice. that the equivalent Q-factor Q™ for the rotor is negative when the rotor 
speed jd is greater than the whirling speeded. Physically this means that the 
damping force F. acts in the direction of whirling for subsynchronous vibration 
and is thus destabilizing. 


Incorporating the expression for the rotor Q-factor into equation A6 yields the 
following equation for the modal Q-factor which includes the effects of spinning 
and whirl ing rotors. 
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Q » 


na PE, nb^ 

n - a -+iT(i-v«) 

a=l Q, b=l D 



(A10) 


PE 

na 

PE a 

nb 
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PE b 
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total system modal potential energy 

number of lumped damping components (mounts, bearings, dampers) 

modal strain energy for lumped damping components 

component Q-factors for lumped damping components 

number of field damping components (casings, frames, rotors) 

spin speeds for field damping components (0, is zero for static 
components) 

modal strain energy for field damping components 
component Q-factors for field damping components 
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